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SEMICLASSICAL APPROACH TO CAPTURE COLLISIONS 
BETWEEN IONS AND POLAR MOLECULES 1 
by John V. Dugan, Jr. and John L. Magee 2 
Lewis Research Center 

SUMMARY 

Quantum mechanical average cross sections are derived for low velocity capture 
collisions between ions and polar molecules. These cross section expressions incorpor- 
ate the rotational energy of the molecule and are derived from the Stark interaction poten- 
tial. 

A classical upper limit which is good for all ion velocities is presented for the cap- 
ture cross section. Also, an approximate cross section is derived from classical con- 
siderations for the regime of high ion velocity by simple averaging over the static orien- 
tations of the polar molecule. 

Finally, a numerical calculation of the capture cross section is performed by random- 
izing initial conditions for two hypothetical molecules with fixed rotational energy at high 
values of ion velocity. The Lagrangian equations of motion are solved to obtain the cap- 
ture cross section from a collision ratio. These cross sections are decreasing functions 
of molecular rotational energy and ion-molecule relative energy. The cross sections are 
discussed in the context of experimental results and the quantum theoretical predictions. 

This study of collision between hypothetical molecular species is intended to serve 
as a preliminary description of such events. 


*A brief summary of some of the results in this report was given at the 17th Annual 
Gaseous Electronics Conference, Atlantic City, New Jersey, Oct. 14-16, 1964 and at 
the 4th International Conference on the Physics of Electronic Collisions at Quebec City, 
Canada, Aug. 2-6, 1965. More extensive details concerning aspects of this report are 
included in M A Semi- Classical Theory of Capture Collisions between Ions and Polar 
Molecules, M John V. Dugan, Jr., Ph. D. Thesis, Dept, of Chemistry, Univ. of Notre 
Dame, Notre Dame, Ind. , May 1, 1965. 
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Professor of Chemistry and Associate Director of Radiation Laboratory, Dept, of 
Chemistry, Univ. of Notre Dame, Notre Dame, Ind. 



INTRODUCTION 


The attractive ion- molecule (IM) collision has been extensively studied because of its 
proposed role in the reaction mechanisms of irradiated gases (ref. 1) . Most cross sec- 
tion data for specific IM collisions have come from mass spectrometric experiments. 
These investigations differ from conventional crossed beam experiments in that average 
cross sections are measured rather than monoenergetic quantities. The observed cross 
section Q is a function of the ion energy e ^ which is proportional to a distance <L 
through the target gas (ref. 2) . The capture cross section is not measured directly but 
has a lower limit set by the value of the reaction cross section (ref. 3) . 

Until recently, the theory employed to calculate IM capture cross sections was iden- 
tical to that used to describe ion-atom collisions (ref. 4). The Langevin polarization po- 
tential 
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had been used to predict the cross section. This potential term describes the interaction 
between an ion and the dipole which it induces in the electronic cloud of the molecule 
(ref. 5). In equation (1), a is the average electronic polarizability of the molecule, 
e the electronic charge, and R the distance between ion and molecule. (A complete 
list of symbols is included in appendix A. ) The cross section for capture was calculated 
from the equations for conservation of energy and angular momentum. For an IM rela- 
tive energy e (at infinite distance), the Langevin cross section is (ref. 5) 


a L 



( 2 ) 


The Langevin theory has proven inadequate in predicting the capture cross section 
for collision between an ion and a molecule with a permanent dipole moment. An addi- 
tional IM term accounting for the permanent dipole interaction is (ref. 6) 
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-pe cos 9 



( 3 ) 


where ju is the dipole moment and 9 is the angle between the negative half of the dipole 
and the R vector. 


Theard and Hamill (ref. 6) included the Vp term with cos 0 = 1 (for simplicity) and 
derived the following cross section contribution: 



( 4 ) 


They then explained the large cross sections observed for certain ion- polar molecule 
collisions in terms of a total a given by + c r^. These experimental results have 
been extended by Moran and Hamill (ref. 7). Their experiments have demonstrated that 
the ion-dipole term does make a large contribution to the observed cross section. The 
choice of equation (4) as the ion-dipole term presumes the dipole can always adjust to a 
maximum orientation in the ion field. This, however, is not generally the case. The ion 
transit time and rotational period would be expected to be important considerations. 

The main purpose of this report is to present a preliminary description of a semi- 
classical treatment of ion-polar molecule capture collisions. It is intended that this ap- 
proach be particularly useful to the experimenter. Throughout the study, the interaction 
is presumed to take place between an ion with velocity v^ and a polar rotator which has 
no translational motion. The theoretical treatment is developed for both linear and sym- 
metric top molecules and is comprised of the following three parts: 

(1) Capture cross sections for low ion velocity are derived from quantum mechanical 
interaction potentials by simple averaging over rotational quantum numbers. In this re- 
gion the ion field is effectively static and the EM interaction is treated as an M ion- 
induced" Stark effect. The manner in which various rotational distribution functions can 
be incorporated to calculate average cross sections is illustrated. 

(2) Classically derived expressions for the capture cross section are presented for 
high ion velocity. This is a region of particular interest for the experimenter. 

(3) Numerical results for the cross section are presented for hypothetical molecules 
of assigned dipole moments and electronic polarizabilities. These calculations are per- 
formed for randomized initial conditions at high ion velocity. 

Finally, the numerical results are compared with both the low velocity predictions 
and the high velocity limiting values. 


QUANTUM THEORETICAL CROSS SECTION AT LOW ION VELOCITY 
Region of Theoretical Validity 

Static field assumption . - The polar molecule Stark effect is defined as the splitting 
of rotational energy levels in a static electric field (ref. 8). In order for the ion field to 


3 



be considered effectively static, its rate of change must remain small compared to the 
molecular rotational frequency. This condition can be written (ref. 9) as 

— — « 1 (5) 

F dt 

o 

In condition (5), F is the ion field (=e/R ) and t R is the rotational period of the polar 
molecule. 

Condition (5) indicates the limit of validity for the low velocity treatment. The nor- 
malized rate of change of the radial component of the ion field is 

JL ^ .1® (6) 

F r dt R 

If the ion velocity is approximately equal to the radial velocity, then condition (5) can be 
written as 



*R <<: 1 


( 7 ) 


-40 

At a rotational energy E D of 0. 05 electron volt, for a moment of inertia of I = 10 

o Jtv 

gram- centimeter (e.g. , hydrochloric acid molecule, HC1), the rotational period is 


t 


R 



1. 57x10' 13 sec 


( 8 ) 


If a 10-percent change in the ion field is allowed during the time t R , the maximum ion 
velocity that satisfies condition (7) is 


w= 318x1 ° 11r c 


(9) 


where R c is the distance of closest approach, which can be chosen as roughly equal to 
the impact parameter b c> This is a lower limit choice for R c , but the choice is not 
important since condition (5) is not the most restrictive condition for the validity of the 
treatment. For R £ = 15 angstroms, the maximum velocity from equation (9) is 4. 8x10^ 
centimeters per second. This value of R is estimated self- consistently by assuming 
ttR = ov» + ov for representative values of dipole moment and reduced mass. This ve- 

C JJ Lt 
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locity corresponds to several times thermal energy for an IM pair with a relative mass 
of 50 atomic mass units. In general, condition (7) will not be fulfilled in the range of ex- 
perimental ion energies; however, it will prove valuable to compare the trends of this 
theory with numerical results. 

Pertur bation approximation . - An additional restriction on the quantum mechanical 
approach arises from the Stark effect approximation. The derivations of the Stark inter- 
action energies between ions and polar molecules are valid only if the field energy is a 
perturbation on the rotational motion (ref. 8). In fact, the ion-dipole interaction energy 
must be much smaller than the separations between rotational energy levels. This crite- 
rion can be written as 


e r 


( 10 ) 


In this expression, J is the rotational quantum number and is the rotational energy. 
The average energy level spacing is about E^/J. For representative values, J = 10, 
ju = 1. 5 Debye units, and Ej^ = 0. 05 electron volt, condition (10) with JpF/E^ = 0. 1 is 
satisfied for an electric field F corresponding to R = 67 angstroms. This IM separa- 
tion is several times the R„ employed previously. In fact, the ratio of the second to 

c o 1/4 

the first R c value can be written as (2Mjiie/h ) ' , which is always greater than unity. 
Thus condition (10), which specifies a minimum R, is more restrictive than equation (9). 


Critical Impact Parameter 


Gioumousis and Stevenson (ref. 4) have carried over Langevin's study of ion trajec- 
tories to calculate reaction cross sections. The interaction potential assumed is equa- 
tion (1) but a critical reaction radius r c is proposed for IM pairs. In reference 4 it is 
shown that ions which start toward the molecule at an impact parameter less than a criti- 
cal impact parameter b c will be captured. These captured ions will spiral inward until 

ion and molecule react inside a critical radius r . This quantity r should not be con- 

c c 

fused with the R c of this report. The latter is the distance of closest approach for the 
impact parameter b c - The critical impact parameter for the Langevin case is 


b 


c 




( 11 ) 


In this study the low velocity trajectories will be treated as in reference 4 for the 
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CjR + C£R potential. An ion which approaches the polarizable polar molecule ex- 
actly at b will enter an unstable circular orbit of radius R . The derivation of the 
c c 

capture cross section from the assumed interaction potential will now be described. 


Interaction Potential 

Quantum mechanical expressions of the rotational Stark effects will be employed to 
obtain terms for the interaction potential. The quantum mechanical derivations of the 
first- and second-order Stark effects appear in reference 8. The total interaction poten- 
tial is 


V int = aW l + iW 2 + V L < 12 > 

where AWj and AWg are the first- and second-order Stark energies, respectively. 

The Langevin term Vg is included because polar molecules have polarizable distribu- 
tions of electronic charge. 

First-order Stark effect . - The interaction energy for the first-order Stark effect for 
the symmetric top molecule is (ref. 8) 


AW 


IT " 


-li FMjK 
J(J + 1) 


(13) 


For the linear molecule, this effect is zero. Equation (13) gives the energy change in a 
static field for a symmetric top molecule such as methyl cyanide, CH^CN. In equation 
(13), J is the total angular momentum quantum number, Mj denotes the component of J 
in the field direction, and K is the projection of J on the symmetry axis of the top. If 
the molecule can orient randomly relative to the field, then Mj ranges from +J to -J 
with equal probability. The quantum number K can also take on all values from +J to 
- J but not with equal probability. The relative probabilities of various K values depend 
on the molecular geometry and gas temperature (ref. 8) . To predict these values one re- 
quires a relation between both J and K and the rotational energy. 

The rotational energy of a symmetric top molecule is (ref. 10) 


E 


RT 


_ J(J+ 1)K 2 , K 2 /h 2 


21 


B 


21 , 



(14) 


In equation (14), Ig is the moment of inertia about an axis perpendicular to the symmetry 
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axis and I A is the moment of inertia about the symmetry axis. Since (for constant tem- 
perature) the energy Eg,p depends on J and K, once J is specified the values of K 
are governed by the molecular geometry. For a prolate top ("cigar like”), I A < Ig, 

E Drr increases with K for constant J so that small K values are most probable. For 
an oblate top ("pancake like’ T ), I B < I A , E RT decreases with increasing K so that large 
values of K are most probable. The distribution of K values is further complicated by 
weighting factors that arise from nuclear spin statistics (ref. 8). 

Classically, the first-order Stark effect is the orientation of a polar molecule such 
that it has a net dipole moment in the field direction. This occurs for a symmetric top 
molecule when the dipole moment has a component along the total angular momentum 

vector (m k/V j ( j + 1)). The quantity Mj, equal to Vj(J + 1) cos 0, is the projection of 

J along the field and the projection of i± is juMjK/J(J + 1). 

For a linear polar molecule (e. g. , HC1), the first-order Stark effect vanishes since 
K = 0 in equation (13). In this case there cannot be a component of j u parallel to J 

since k/Vj(J + 1) is zero. This is strictly true only for diatomic molecules. It does 
not apply to linear polyatomic molecules which are excited to bending modes (ref. 8). 

The rotational energy of a diatomic molecule is (ref. 11) 

E (15) 

R 21 


where I is the moment of inertia about an axis perpendicular to the internuclear (polar) 
axis. 

Second- order S tark effect . - The second- order Stark interaction energy for the sym- 
metric top molecule is (ref. 8) 


AW 


2T 


2„2 T 

p F I B 


h 


(J 2 - K 2 )(j 2 - Mj) [(J + l) 2 - K 2 ] [(J + l) 2 - M 2 


(16) 


«r(2j - i) (2j + i) 


(J + ir(2J + 1)(2J + 3) 


For the linear molecule, K = 0 and AWg becomes 


o o r ox 

I R p F J(J + 1) - 3M^ 

AW 2R = < 

H 2 J(J + l)l< 2j - 1 )< 2j + S )J 


(17) 


Classically, this effect corresponds to a nonuniform rotation of the molecule in the 
electric field. Under such conditions the dipole may spend most of the time lined up 
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either with or against the field. Thus the molecule will exhibit a net attraction or repul- 
sion for the ion. Quantum mechanically, this corresponds to polarization of the wave 
function of a rigid rotator (ref. 8) . 


Capture Cross Section for a Linear Polar Molecule 

Case of constant rotational energy . - The potential for the interaction between an ion 
and a linear polar molecule is 

( V lnt) R = iW 2R +V L 


1 

-2/i 2 I 

r j(j + 1) - 3m 2> 

2 a e 2 
e o 

2 

fi 2 J(J + 1) 

[(2J - 1)(2J + 3)^ 

R 4 2R 4 


Since the square of the electric field appears in AW 
interpreted as a rotational polarizability term a iR . 
is defined as 


2 R , the quantity in brackets can be 
If then the total polarizability a^ R 


“tR * “iR 


+ ot. 


equation (18a) can be written as 


a 


tR " 



(18b) 


(18c) 


Equation (18c) is analogous to the Langevin expression of equation (1). The partial cross 
section for a given Mj value can thus be written in the same form as the Langevin cross 
section. The electronic polarizability a Q is merely replaced by the total polarizability. 
Positive polarizabilities correspond to attractive interactions; repulsive collisions do 
not contribute to the capture cross section. The partial cross section is 


°tR ~ ffP J 


2c W M j) e ' 


1/2 


(19) 
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Relative number of molecules, N .( J) 




(b) Symmetric top polar molecule (methyl chloride) at temperature of 
300° K (ref. 10). 


(a) Linear polar molecule (hydrochloric acid) at temperature of 300° K 
(ref. 11). 


where P T is the probability that the mole- 
cule is in the ±Mj state. The term 
a tR * s a f 11110 ^ 011 °f Mj, the projection 
of J along the ion field. In equation (19), 
e = Mvj/2 is the relative translational 
energy of the IM pair written in terms of 
the ion velocity v.. The reduced mass 
M is M q Mj/(M o + Mj) where M Q and 
Mj are the molecule and ion masses, re- 
spectively. 

Molecules with small moments of in- 
ertia will be in relatively low J states at 
normal gas temperatures. For example, 
the distribution of molecules in J states 
for HC1 at 300° K is shown in figure 1(a) 
(ref. 11). For such "quantum rotators, M 
the cross section must be calculated by 
summing all positive partial cross sec- 
tions. Their values range from the first 
integral value of Mj which gives +a^ R 
up to J. 

Molecules in high J states, that is 
"classical rotators, " have distributions 
as shown in figure 1(b). For these mole- 
cules, the M t states are effectively con- 


Figure 1. - Distribution function among states of rotational quan 
turn number J. 


by integrals. 


tinuous and the summations over the posi- 
tive polarizability values can be replaced 
As shown in appendix B, in neither case does the cross section distinguish 

This is because the second-order Stark effect does not 

8 ). 


between positive or negative Mj. 
remove the M T degeneracy (ref. 


The average cross section for zero electronic polarizability at high J taken from 
appendix B is 


, 2 2 
— „ P e 

C-TD = 7T l — 

lR '4E r € 


1/2 


V6 - ln(V3 + Vi) 

Ve 


= 7T 



2 2 


1/2 


(0. 53) si a 


2 l 4E R e 


( 20 ) 
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IM relative energy, € , eV 


Figure 2. - Variation of average quantum mechanical cross section with relative ion-molecuie energy for 
various molecular rotational energies. Linear polar molecule with zero electronic polarizability; 
dipole moment, p = 1. 67 Debye units; calculated from equation (20). 


This cross section is plotted against e in figure 2 for five different rotational energies 

and a constant value of p (a Q = 0). The upper limit to the region of exactness of the 

perturbation approximation is indicated by the solid curves. In all subsequent figures, 

dashed curves are used to indicate extensions beyond the ranges of validity. 

The average cross section for nonzero electronic polarizability is where 
2 

a Q > p /4Ej^ from appendix B: 


6 \ .rT . — ; 
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IM relative energy, e, eV 

Figure 3. - Variation of average quantum mechanical cross section with relative ion-molecule energy 
for various molecular rotational energies. Linear polar molecule with nonzero electronic polariz- 
ability, a Q « 1.67 angstroms 3 ; dipole moment, p = 1. 67 Debye units. 


For Q' o <fx 2 /4E R , 



( 22 ) 


The cross sections from equations (21) and (22) are plotted against e in figure 3 for 

Q 

jut = 1. 67 Debye units and a = 1. 67 angstroms . 

It must be emphasized that these average cross sections have been derived for fixed 
values of rotational energy (i. e. , constant J). Such a procedure ignores the distribution 
of J values at the equilibrium gas temperature. These expressions do, however, prove 
valuable for comparison with numerical results even though the curves of figures 2 and 3 
are extended considerably beyond the range of validity. 

Rotational energy distribution funct ion - linear molecule . - In equations (20) to (22) 
representative rotational energies were chosen that would be equal to kT for five values 
of temperature from 116° to 11, 600° K. These values of range from 0. 01 to 1. 0 
electron volt. 
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In an experiment, however, the target molecules are actually distributed over a 
range of rotational energies as described by the distribution function f(J) . This normal- 
ized function is (ref. 11 ) 


f(j) = (2J + !) e - J ( J+1 ) Bh /kT ( 23 ) 

kT 

The average cross sections at a temperature T for constant e are the normalized inte- 
grals of equations (20) to (22) weighted by f(J). 

Capture Cross Section for Symmetric Top Polar Molecule 

Case of constant rotat ional energy. - For the symmetric top, the first-order Stark 
term (eq. (13)) must be included in the total interaction potential, in addition to 

the second-order Stark term AW 2 rp (from eq. (16)). The electrostatic interaction poten- 
tial is then 



-jueMjK 

R 2 J(J + 1) 



(24) 


where 


2 AW 9t 

V = + 


is the total polarizability of the polar top. This quantity is the sum of the induced and 
electronic polarizability terms (R~^ terms). The partial cross section for capture can 
be derived by following an approach outlined in reference 5. 

The effective potential for the ion- symmetric top interaction is 

2 

v efT + ( v tat) T < 25 > 

2MR^ 1 

2 2 

In equation (25) the centripetal contribution L /2MR is a fictitious potential term where 
L is the angular momentum (Mv^b). This term accounts for the relative motion of the 
IM pair, and its inclusion in this manner allows for a straightforward cross section 
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derivation. In conventional notation, the first-order interaction term AW.^ may be 
written as -ju^e/R 2 where / u 12 is the dipole moment matrix element (ref. 8). Let the 
effective angular momentum term be 



L 2 - 2n 12 eM = (Mvjb) 2 - 2]u 12 eM 


(26) 


The maximum of the effective potential is at the critical distance of closest approach on 
the crest of the angular momentum barrier. Since 


then 



(27a) 


T 2 2 
L eff = a tT e 

2 MR 2 R 4 
c c 


(27b) 


From equation (27b), the expression for the distance of closest approach is 


R 


2 

c 


2Ma tT e 


J eff 


(28) 


The translational energy at infinite distance is equal to the effective potential at R . 
This is true if the rotational energy of the molecule is effectively constant. This condi- 
tion can be written as 


e 



However, 



is simply 


1/2(L 2 ^/2MR 2 ^ because 


eff 


ef f 
2 MR 2 


a tT e _ L eff 

2R 4 2MR 2 
c c 



(29) 


(30) 
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The final form of equation (30) is obtained by substitution for c^e /R from equa- 
tion (27b). 

Substituting the expression for R 2 from equation (28) into equation (30) gives the 
following relation for e : 

L eff _ - 2 >*12 eM ] 2 (3 

8M 2 a tT e 2 8M 2 a tT e 2 

In equation (31), b c is the critical impact parameter corresponding to the distance of 
closest approach R c> Rearranging and taking the square root of both sides of equation 
(31) give 


(Mvjb c ) 2 - 2ju 12 eM = (8M 2 a tT e 2 e) 


1/2 


(32) 


If equation (29) is used for Mv 2 /2, the partial cross section for the ion-symmetric top 
capture collision is 


°tT “ * h c ~ 57 




''lOL.mQ 

tT 


, 1/2 


(33) 


In equation (33) both terms are functions of Mj and K. 

The partial cross sections for the symmet- 
k ric top cannot be simply averaged over the 

quantum numbers Mj and K. Although all 
values of Mj are equally probable, the values 
of K for constant J obey a distribution func- 
tion fj(K). However, it is more conventional 
to consider the distribution function f^.(J) for 
constant K. 

Distributions f K (J) for the parameter K 

are shown in figure 4 for a sample symmetric 

top. For this molecule it is assumed that all 

nuclei have very large spins. Levels for which 

K is a multiple of three must be doubly weighted 

Figure 4. - Distribution function among states of rotational f or particles of spin 1/2 which rotate about the 
quantum number J for a symmetric top polar molecule 

for various quantum numbers, K, and l B • 5I A symmetry axis. The K value assignment to be 

(ref. 10). 
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made in the next paragraph was consistent with the initial conditions of the numerical 
study. 

The partial cross section for the ion- symmetric top capture collision at high J for 
a fixed K is obtained by substituting for an< * a tT equation (33). This partial 

cross section is 


2 2^/ 2 

tT ^ J 4J \ E R e 


2\/2Je 


J 2 ^ . 5 ] + M 2 


19 


1/2" 


(34a) 


In equation (34a) the value of K = ±j/2\j2 is assigned by solving equation (14) with 
I R = 51^. This moment of inertia ratio was used in the numerical study. 

The average cross section is 


°tT - 


7 Tjie 

2V2J 2 < 


ax Mj J ^ax ^2 ± & 2\ 1/2 , 

^(Mj) ^\ Er€ J J(m,\ 


mm 


K) . 


mm 


(34b) 


In equation (34b) 


16E 


R 


- 5 


2 e c 


a = 


19 


i.1* 

and the probability that the top is in the +Mj state Pj = (1/2 J) has been included. 
Equation (34b) integrates to 


J tT 


_ 7rpe 


°* nri y- 2 \ max M min ) 


2 \ 
miny 


\eJ 


/1Q 2 2 X 1 / 2 r 

n j 19/j. e 


8J 


E R e 


| M j( 3 


m 2 2N 1 / 2 2 

M J ± a ; • - 


± a In 


[ m j + (> 


2 . _2\ 1 / 2 I 


M j ±a V 


( m j). 

( m j) 


max 


min 


(34c) 
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Average cross section, o^j, A 2 


0 2 4 6 8 10 

Quantum number, K 

Figure 5. - Distribution function among states 
of quantum number K for symmetric top 
polar molecule from figure 4 for rotational 
quantum number J of 9. 



IM relative energy, € , eV 


Figure 6. - Variation of average quantum mechanical cross section with relative ion- 
molecule energy for various molecular rotational energies. Symmetric top polar 
molecule with nonzero electronic polarizability, a 0 = 1.67 angstroms 3 ; dipole 
moment, |i = 1. 67 Debye units; calculated from equation (34c). 
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The limits of integration in equation (34c) are defined by the condition that the integrand 
of equation (34a) be positive. 

The rotational energy is again held constant and corresponds to a representative 
molecule with Eg,p = (3/2)kT. For simplification it is assumed that the initial K value 
is specified by equipartition of energy and is ±<1/2^2. Actually the distribution fj(K) for 
constant depends on the molecular geometry and temperature. A sample of the 

function fj(K) for a top with Ig = 51^ is shown in figure 5. The average cross section 
for fixed rotational energy and nonzero electronic polarizability is plotted in figure 6. 

Rotational energy distribution function - sy mmetric top . - The average cross section 
for constant J and a fixed K value is not directly usable by the experimentalist. Rather, 
the average over the total distribution function f(J, K) is the experimental quantity of in- 
terest. This normalized distribution function is (ref. 8) 


f(J,K) 


B 2 Ah 3 


[_7r(kT)' 


1/2 


(2j + 1)e -[J(J+l)B+K 2 (A-B)]K 2 /kT 


where A = h/W 2 I^. As mentioned earlier, the moment of inertia 1^ will be greater 
than Ig for an oblate top and less than Ig for a prolate top. For a constant tempera- 
ture the lower values of K are more probable for a prolate top. The higher values of K 
become more probable for an oblate top. 

The average cross sections at a given e which are of interest to the experimentalist 

are 


u(e) = f f a(J, K)f(J, K)dJ dK 


(35) 


The cross section a(J, K) must be previously integrated over a random distribution of 
Mj values. 


CLASSICAL CROSS SECTIONS AT HIGH ION VELOCITY 
Two Limiting Cross Section Expressions 

In the region of high ion velocity, two classical expressions can be obtained for the 
capture cross section for extreme conditions of collision behavior. The first is a rigor- 
ous upper limit to the capture cross section obtained from the condition that the dipole 
can always orient so as to attract the ion. The second is an approximate cross section 
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Figure 7. - Variation of classical maximum cross section with relative ion-molecule 
energy. Nonzero electronic polarizability, a Q = 1.67 angstroms^ dipole moment, 

(j = 1. 67 Debye units; calculated from equation (36). 


obtained from the condition that the collision time is much shorter than the rotational pe- 
riod. The latter expression is not rigorous but is conceptually useful in discussing colli- 
sion behavior in this regime. 

Upper limit cross sec tion. - As explained earlier, the maximum cross section for 
an ion-polar molecule collision independent of molecular geometry is 


a 

max 


TT/xe 

€ 


+ cr T 


(36) 


This is precisely the expression used in references 6 and 7 to fit experimental IM data. 

It says that regardless of the ion trajectory the dipole can always orient favorably for 
capture. This is a type of adiabatic approximation where it is assumed that the Hamil- 
tonian changes very slowly. This rate of change must be sufficiently slow so that the dy- 
namical system can always act to minimize the energy. Such an approach clearly ignores 
the rotational energy of the target molecules. This assumption actually violates conser- 
vation of energy since it takes no account of rotational acceleration or deceleration ac- 
companying collision. This cross section is plotted in figure 7. 
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Estimated capture cross section . - The condition for the static orientation approxima- 
tion is that the IM interaction time be very much shorter than the period of molecular ro- 
tation (ref. 12). This condition can be written as 



(37) 


where t is the collision time and t„ is given by equation (8). The quantity t can 
best be written in terms of the relative velocity and a characteristic interaction distance 
R^. A meaningful choice for R^ is the distance at which the rotator can first be com- 
pletely ’'hindered” by a static ion field, that is, cannot make a complete rotation. Then 
R h is defined by 


R h = 



(38) 


The quantity t becomes 


-14 

For a molecule with = 0. 05 electron volt = 8x10 erg, the interaction distance for 
a dipole moment of 1. 67 Debye units is 10 angstroms. The minimum velocity at which 
the sudden approximation is valid, for tjj/t c = 10 and t^ given by equation (8), is 



min 


10" 1 R, 


R 


10 


-8 


1. 57X10 


•13 


— 1 6 — 1 
cm sec” = 6x10 cm sec 


(40) 


Letting M cos ® and = a Q in equation (33) yields the partial cross section 

Oq = 7 rb^ = cos 6 + ct l (41) 

Equation (41) is written with the assumption that the form of equation (33) is valid both 
for the high and low interaction energies (velocities). Equation (41) means that all ions 
which approach at less than the critical impact parameter b c will be captured. 
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It can be argued that some ions which are initially attracted as specified by equation 
(41) will enter orbits that later become repulsive. This effect could occur as the ion 
spirals inward but before the R value becomes less than the r c of reference 4. (The 
numerical results indicate only a small fraction of collisions behave this way. ) The 
averaging procedure which will be described does, however, provide an instructive limit- 
ing case. 

The average cross section is simply the normalized integral over the partial cross 
sections independent of rotational energy. It is assumed that the dipole has no time to 
adjust in the ion field. The differential solid angle is df2 = sin 9 d9 dcp . The static ori- 
entation average cross section is 



7 Tj ne 
e 



e 

max 

cos 9 sin 6 dd 


(la e 2 \ /*Vax 
+ £ f — 2 — j j sin 9 d 9 


(42a) 


s * 9tt 

The factor of 1/2 in equation (42a) comes from the / " d cp divided by the solid angle 

normalization factor 4m . The value of 0 is determined from the condition that 

max 

equation (41) be positive. In terms of cos 9 this occurs when 


< cos 0 )min = - 


'2 a e 




1/2 


, 1/2 


(42b) 


where the critical energy e c is (y?/2a^j. 
Equation (42a) can be written as 


- True 
6 2e 


»/(cos 9) 


(2 a e 

cos 9 d(cos 9) + — ' ° 

2 


2 \l/2 


mm 


/ 


d(cos 9) 
(cos 9) 


(43) 


'min 


Upon integration, and substitution of the expression for (cos 0) m ^ n » equation (43) becomes 


a 

6 4e 



(44) 


20 




1 CT 4 10‘ 3 10' 2 10' 1 10° 


IM relative energy, e, eV 

Figure 8. - Variation of average static orientation cross section with relative ion- 
molecule energy. Nonzero electronic polarizability, a 0 = 1.67 angstroms^; 
dipole moment, \i = 1. 67 Debye units; calculated from equation (44). 

This cross section of equation (44) reduces to cr-^ = 77^2o! o e V ^ cp itical energy 

e . It is plotted against e in figure 8 for jjl = 1. 67 Debye units and a = 1. 67 ang- 
3 o 

stroms independent of rotational energy. Condition (37) is, of course, the lower limit 

to the region of validity of this sudden approximation. From equation (40) for an E R 

value of 0. 05 electron volt and M = 50 atomic mass units, the minimum value of e is 

1000 electron volts. This minimum energy is three orders of magnitude greater than the 

largest value of the abscissa in figure 8. 

Numerical Treatment 

Need for a numeri cal high velocity treatment . - Quantum rotators have energy level 
separations which are significant compared to kT and cannot be considered continuous. 
However, the IM interaction potential at high ion velocity cannot be considered a pertur- 
bation on the rotational motion. In fact, neither J, K, nor Mj are good quantum num- 
bers when critical impact parameters are of the order of molecular diameters. The 
complete approach to a collision theory in this regime, therefore, must be by means of 
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a quantum mechanical treatment. Such a description would require an exact calculation 

of a time -dependent strong field Stark effect. That problem is clearly intractable and its 

consideration is beyond the scope of this report. On the other hand, classical rotators 

with J about 10 or greater have energy levels that can be regarded as continuous. Most 

-40 2 

molecules with I values greater than 10 gram- centimeter will be excited to these J 
values at normal temperatures. From this viewpoint the high velocity IM interaction can 
be considered purely classical. However, both the direction and magnitude of the rota- 
tional angular momentum will vary during interaction and these effects should also be 
considered. To treat this classical collision satisfactorily it is necessary to solve the 
equations of motion numerically for the entire IM system. This method of attack will be 
outlined in the next section and results will be presented for hypothetical molecules. 

Models and general approach . - The ion- dipole interaction in the classical limit was 
studied by numerical solution of the equations of motion for two different models, which 
are approximations to real molecular systems. The target for model I was a rigid rota- 
tor polar molecule with a permanent dipole moment along the internuclear axis. This 
polar molecule was allowed to interact with an ion moving translationally relative to its 
center of mass. The target for model II consisted of a rigid symmetric top molecule 
with a dipole "frozen" in it along the symmetry axis. The ion was treated just as in 
model I. 

Coupled differential equations were written for the two physical models, and the 
models were programed in Fortran IV for solution on the IBM 7094 digital computer. 

The Lagrangian formulation of the force equations was utilized for each model. The con- 
servation of energy equation was carried along to check the accuracy of the integration 
technique. 

The coordinate system used for model I appears in figure 9. A discussion of this 
system is included in appendix C. The values of the constants chosen for the linear mole- 
cule and a discussion of the translational energy approximation appear in appendix D; 
also included in appendix D is a discussion of the assignment of initial conditions to the 
variables and their time derivatives. 

A total of 36 cases was solved for constant (initial) rotational energy for each impact 
parameter. The orientations of the negative end of the dipole and plane of rotation were 
varied relative to the ion path. This was done to completely randomize the initial condi- 
tions on the IM pair with the ion effectively at infinity. 

The IM interaction potential is composed of both the ion- dipole and polarization 
terms. The polarization term is simply -« o e /2R . However, the ion-dipole term is 
a complicated expression containing four different angles. The IM total energy expres- 
sion and derivation of the equations of motion are also included in appendix E. 

The numerical study of the symmetric top target molecule was performed to investi- 
gate the effect of molecular geometry on the time- dependent IM interaction. The coordi- 
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Figure 10. - Coordinate system employed in computer study of interaction between an ion and a symmetric 
top polar molecule. 



nate system used to describe the ion-symmetric top interaction is shown in figure 10. 
The randomization of initial conditions was performed similar to model I. Additional 
considerations regarding the extra degree of rotational freedom are discussed in appen- 
dix F. The derivation of the equations of motion is also included therein. 

Definition of the numeri cally calculated cross section. - The numerically calculated 
capture cross section is 


ci 



(45) 


where i = 1 or 2 identifies the molecular model. The critical impact parameter b c - 
which defines the capture cross section is not a conventional impact parameter like b c . 
The latter is the value of b below which capture collisions occur with unit probability. 
However, although the -a o e /2R potential is always attractive, the ~/j.e/R cos 8 
term has a repulsive phase. Since the ion permanent dipole potential can thus prevent 
capture of ions, b . is best defined in terms of a collision ratio C D . The quantity C D 
is the fraction of IM cases tried which result in collision for an assigned impact param- 
eter b. The capture cross section is then defined by 


"ci - /° rO> < 46 > 

o 

The quantity dA c in the differential collision area or cross section d( 7 rb ) = 2nb db. 
Equation (46) can then be rewritten as 


ci 


= 7 T 



C R (b) d(b 2 ) 


(47) 


where bp is the largest value of b for which any collisions are observed. The square 
of the critical impact parameter for capture is simply 

b ci = r° C R (b) d(b2 ) ( 48 ) 

Jo 

/ 2 2 
b db = b /2 and the 2 cancels the 2 from d(b ). 

The square of the critical impact parameter for capture is obtained directly by measure- 

2 

ment of the area under the curve C R against b . 
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RESULTS AND DISCUSSION 


Qualifications Regarding the Quantum Theoretical Cross Section 

-38 

For the hypothetical molecule examined numerically, 1 = 3. 85x10” gram- 
2 

centimeter and the IM interaction cannot be considered a perturbation for velocities 
greater than 10 centimeters per second. This value of Vj corresponds to an IM energy 
of about 10” ® electron volt for a representative value of reduced mass. This velocity is 
1/50 of the minimum value assigned in the numerical calculation. In the regime of veloc- 
ities included in the computer study, transitions between rotational states do not obey 
quantum mechanical selection rules; rather, the rotational energy varies through a con- 
tinuum of states under the influence of the ion field. This field is time dependent and a 
quantum mechanical description of the IM interaction is an intractable problem. However, 
it does prove useful to assume that the continuum average of the matrix elements approxi- 
mates the classical time-averaged behavior. There is no guarantee that this situation will 
hold true in the region where numerical results are to be compared with theoretical pre- 
dictions. Nevertheless, the average quantum mechanical expressions do provide a con- 
venient framework for qualitative discussion of numerical results. 

It is also true that the quantum mechanical cross sections were derived for constant 
J. It is obvious that J and Mj will vary during a representative interaction. In the 
derivation of the quantum mechanical expressions it was assumed that all values of Mj 
remain equally probable for a fixed J. This approximation is equivalent to assuming 
that Mj varies randomly during an average collision. 


Numerical Results for Model I 

o 

Curves of C D against b for E D = 0. 043 and 0. 086 electron volt are shown in fig- 

rt JrC 4 5 

ures 11(a) and (b). These plots are for ion velocities of 5. 5x10 and 10 centimeters per 
second, respectively. These velocities correspond to IM relative energies of 0. 05 and 
0. 16 electron volt for an M value of 32 atomic mass units. The shape of such plots is 
sensitive to values of both rotational energy and ion velocity. The error in calculating 
the numerical cross section is roughly b c ^ In, where n denotes the number of cases at- 
tempted for random initial orientations at a given b value. 

Computer checks of the orientation of the rotational plane indicate that it is effec- 
tively random for a representative collision. This gives credibility to the assumption 
that all Mj values were given equal probability in the derivations of the quantum me- 
chanical cross sections. The cross section values a c j against e are shown for the 
two rotational energies together with the quantum mechanical curves in figure 12. On 
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Impact parameter, b , A* 1 

(b) Ion velocity, 10^ centimeters per second. 

Figure 11. - Variation of collision ratio with free particle impact parameter for two different rotational 
energies. Nonzero electronic polarizability, a n = 1.67 angstroms^; dipole moment, pi = 1. 67 Debye units 



Figure 12. - Comparison of numerically calculated cross section a cl and average 
quantum mechanical cross section 6^ for linear molecule. Nonzero elec* 
tronic polarizability, a 0 = 1.67 angstroms 3 ; dipole moment, p = 1. 67 Debye units. 


the simple comparison basis which has been adopted, the fact that the <x ^ values are 
20 to 80 percent larger than the quantum theoretical values can be accounted for by the 
net lowering of the initial Ej^ value during collision. Monitoring of the machine problem 
was done with an additional subroutine. Results from this check indicate that the rotor is 
usually decelerated when the ion reaches its turning point. The average value of J cal- 
culated in the numerical study can then be interpreted as an indication of the degree to 
which the rotation is polarized by the field. Classically, this effect corresponds to non- 
uniform motion of the dipole in the ion field. Once the ion reaches the hindering distance 
R^ the field tends to line up the dipole but is allowed only a limited time to do so. The 


- 1/2 


dependence which is 


dependence of o c j on rotational energy is very nearly the 
predicted for a^. 

Comparison is made between the cr , values and predictions of tr (eq. (36)) and 
0q (eq- (44)) in figure 13. Below e = 0. 10 electron volt, the numerical cross section 
shows a more gradual dependence on e than the o' curve. The curve of a Q as a 
function of e is somewhat different from the plots of o ^ as a function of e. 

Comparison between o ^ and the Langevin cross section from equation (1) is made 
in figure 14. Both cr c j plots give considerably larger cross sections than cr-^ through- 
out the range of e . It is evident that reaction cross sections measured for linear polar 
molecules are mainly the result of the ion-permanent dipole interaction. This situation 
will hold true unless / u is very small (e. g. , the carbon monoxide molecule, CO). Al- 
though the Lagrangian contains the sum of two potential terms, their average contribution 
to the capture cross section is not simply the maximum cr^ + o^. Meaningful experimen- 
tal cross sections can be obtained by integrating the numerical results over the distribu- 
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Figure 14. - Comparison of numerically calculated cross section a cl for linear 
molecule and Langevin cross section o L . 







Figure 15. - Comparison of numerically calculated cross sections and average 
quantum mechanical cross section for symmetric top molecule. 


tion function in rotational energy. 


Numerical Results for Model II 

The numerical results for the ion-symmetric top capture cross section are shown in 
figure 15. This capture cross section displays the same dependences on rotational and 
translational energy as the linear molecule cross section. 

The average quantum mechanical cross sections are also plotted in figure 15 for 
Ert = 0. 063 and 0. 123 electron volt. The main contribution to the quantum theoretical 
cross section, equation (34b), is from the polarizability term. This behavior is consis- 
tent with the fact that the numerically calculated symmetric top cross sections are only 
slightly larger than the corresponding linear molecule values of figures 12 to 14. These 
results indicate that the additional degree of rotational freedom for the symmetric top 
has little effect on the capture cross section for this hypothetical molecule. Classically, 
this corresponds to a negligible first-order Stark contribution to the capture cross sec- 
tion. The values of the numerically calculated cross sections are smaller than the aver- 
age quantum mechanical values throughout the range of e. This is in contrast to the be- 
havior of the ion- linear molecule capture cross section. 

Since the component of angular momentum along the symmetry axis will vary during 
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interaction, the ion- symmetric top collision has an additional complication. Not only J 
and Mj but also K will take on values different from the initial value as the ion ap- 
proaches the molecule. Any detailed discussion of the similarity between the two sets of 
numerical results must await further numerical checks and interpretation. 


Comparison of Experimental and Numerical Results 

As shown in figures 13 and 15, the numerical capture cross sections for the hypo- 
thetical molecules are smaller than the maximum IM capture cross section of equation 
(36), especially at low values of e. The latter expression has been used to fit data from 
experiments done at gas temperatures of 500° K (refs. 6 and 7). The numerical results 

also indicate the capture cross section has a gradual dependence on IM relative energy 

- 1/2 - 1 
that is less than e ' . This should be contrasted with the e dependence which has 

been reported for certain IM reaction cross sections (refs. 6 and 7). This disparity may 

suggest that the role of spurious effects accompanying low electric field experiments 

should be investigated. Other experimenters have reported difficulty in specifying the 

ion energy and trajectory at low electric fields (ref. 13). 

CONCLUDING REMARKS 

Up to the present, experimental studies of IM capture cross sections have ignored 
the effect of rotational energy of the target molecule. However, the present numerical 
treatment indicates that the capture cross section has an inverse dependence on for 
the hypothetical molecules treated. This suggests that IM reaction cross section data 
should be analyzed with regard to the effects of rotational energy. Also, ion energy de- 
pendencies of such cross sections should be reviewed in the context of the numerical re- 
sults. More detailed comparisons using typical molecules (rather than hypothetical) are 
required before the description of these collisions is complete. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, November 3, 1965. 
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APPENDIX A 


A 


B 

b 


ci 


C 


'R 


E 


R 


E RT 
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SYMBOLS 


2 

rotational constant, h/87r I., 

-1 A 

sec 

collision area in the numerical 
2 2 

study, 7rb , cm 

2 - 1 

rotational constant, h/8jr I, sec 

classical free particle impact pa- 
rameter, L/Mv-i, cm; 

o — ft 

1A = 10 0 cm 


critical impact parameter for 
IM capture collisions (defines 
cross section), cm 

numerically calculated critical 
impact parameter for capture, 
cm 


constant appearing in ex- 

pression, 4E r & 0 /m 2 = 2E R /e ( 
dimensionless 


collision parameter employed in 
computer study, ratio of number 
of collisions observed to number 
of cases tried, dimensionless 

rotational energy of linear mole- 

- 12 

cule, ergs; 1.6x10 erg=leV 

rotational energy of symmetric 
top polar molecule, ergs 

unit electronic charge, 

4. 8x10”^ esu charge units 

2 

electric field of ion, e/R , V/cm; 
300 V/cm = 1 esu field unit 

radial component of ionic electric 
field, V/cm 


f(J) 


h 


K 

I 


J 

K 


k 

L 

L 1 

Lr 

£ 

i 

M 


rotational distribution function, 
probability that the rotator is in 
interval J to J + dJ at J, di- 
mensionless 

-27 

Planck's constant, 6.63x10 
erg- sec 

-27 

h/27T, 1. 1x10 erg- sec 

2 

moment of inertia, (g)(cm ) 

moment of inertia about symmetry 

2 

axis of top, (g) (cm ) 

rigid rod moment of inertia; also 
larger moment of inertia for 

2 

prolate symmetric top, (g)(cm ) 

rotational quantum number, dimen- 
sionless 

projection of J quantum number 
on symmetry axis of top, dimen- 
sionless 

- 1 fi 

Boltzmann's constant, 1.38x10" 
er g/(°K) (molecule) 

classical translational angular mo- 

2 

mentum, Mv^b, (g)(cm )/sec 

Lagrangian function for model I, 
ergs 

Lagrangian function for model n, 
ergs 

variable length of an ion through 
the target gas, cm 

reduced mass, = M q M^/ 

(M, + M ), g; atomic mass 
unit (amu) = 1. 6x10 ^ g 


31 


M, 


M„ 


M 1 


Ni 


N 


JK 


projection of rotational quantum 
number in direction of ion 
field, dimensionless 

mass of polar molecule, g 


RpjEF reference distance 

(ft nondimensional distance 

T gas temperature of target mole- 

cules, °K 


mass of ion, g 

t i 

kinetic energy term for model I, 

relative number density of linear 


ergs 

rotators that are in state 

T 2 

kinetic energy term for model n, 

with rotational energy, 

Cl 

ergs 

e r = J K /2I 

t 

time, sec 

relative number density of sym- 
metric top rotators that are in 

*C 

effective collision time, R^/Vr 

a particular state with quantum 


sec 

numbers J and K and rota- 


characteristic period of rotation 

tional energy, E RT = 

for polar molecule, sec 

j2 h 2 + k 2 /h 2 _ h 2 \ 

^REF 

reference time 

2i b \ 2i a 2i b/ 

V eff 

effective interaction potential for 


probability that rotator is in 
±M t ^ state, 1/(J + 1), dimen- 



sionless 

V int 


probability that symmetric top 
is in a state specified by a 

V L 


particular Mj and a particu- 

V-r, 


lar K, 1/(2J + 1), dimension- 

P 


less 

V 1 

Q 

o 

observed cross section, cm 

V 9 

R 

distance between centers of mass 
of ion and molecule, cm 

£i 

R c 

distance of closest approach for 
the effective potential treat- 
ment at low ion velocity, cm 

V 

max 

R h 

hindering distance at which max- 
imum potential interaction is 
equal to the rotational energy, 

v min 


(pe/Ej^) 1 / 2 , cm 

V 1 


ion- symmetric top system at 
low ion velocity, ergs 

total interaction potential, ergs 

Langevin polarization potential, 
ergs 

ion-dipole potential term, ergs 

potential energy term for 
model I, ergs 

potential energy term for 
model n, ergs 

upper limit velocity for which 
quantum mechanical average 
cross section is valid, cm/sec 

lowest ion velocity at which sud- 
den approximation is valid, 
cm/sec 

ion velocity, cm/sec 
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AW j 

first-order Stark interaction en- 

e 

angle between ion field F and 


ergy, ergs 


negative end of dipole, rad 

AW 1T 

first-order Stark interaction en- 

M 

permanent dipole moment of 

ergy for symmetric top, ergs 


polar molecule, esu-cm; 

aw 2 

second- order Stark interaction 


1 Debye unit = 10" 0 esu-cm 


energy, ergs 

^ 12 

dipole moment matrix element 

* W 2R 

second- order Stark interaction 

for symmetric top, esu-cm 

energy for rigid rotator, ergs 

k 

Euler angle for rotational motion. 

aw 2t 

second- order Stark interaction 


rad 

o 


energy for symmetric top, 

a 

microscopic cross section, cm 


ergs 

CT ci 

numerically calculated collision 

“iR 

induced polarizability of rigid 

q 

2 2 

cross section, Trb^, crn 

polar rotator, cm° 

^cl 

numerically calculated cross 

«o 

average electronic polarizability 

section for model I, cm^ 


of molecular ground state, 

3 

a c2 

numerically calculated cross 


cm° 

section for model II, cm^ 

“tR 

total polarizability of rigid ro- 

3 

°D 

maximum ion- permanent dipole 


tator, cm 

contribution to capture cross 

01 tT 

total polarizability of symmetric 


section, cm^ 

top molecule, cm^ 

°iR 

partial cross section for induced 

€ 

relative (translational) energy of 

polarizability contribution to 


ion molecule pair, ergs 


ion-rigid rotator collision. 

e c 

critical relative (translational) 
0 / 


cm z 


energy, \i /2a Q , ergs 

°iR 

average cross section for in- 

e ! 

translational energy of ion, ergs 


duced polarizability collision 

X 



of ion with rigid rotator polar 

^1 

azimuthal angle for rotational 


molecule, cm^ 

^2 

motion of rigid rod (model I), 
rad 

azimuthal angle for rotational 

°L 

Langevin cross section, 

ff ( 2a o e2 A) ^ » cm2 

motion of the symmetric top 

°max 

maximum IM capture cross sec- 


(model n), rad 

tion, Oj-j + Oj ■ , cm^ 

0 

polar angle for translational mo- 




tion in models I and II, rad 
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II 


°tR 

°tR 


partial cross section for ion-linear 

o 

molecule capture collision, cm 

average cross section for ion- 
linear molecule capture colli- 
sion, cm 2 


cr tT partial cross section for ion- 
symmetric top collision with 
nonzero electronic polarizabil- 
ity, cm 2 


a tT average total polarizability cross 

section for symmetric top, cm 2 

<Jq partial cross section for static 

orientation cross section at 

2 

high ion velocity, cm 


a e 

t r 

$ 

<p 

*2 


a 


static orientation average cross 

section for high ion velocity, 

2 

cm 

nondimensional time 

azimuthal angle for translational 
motion in models I and n, rad 

azimuthal angle for spatial integra- 
tion, rad 

coordinate for rotation of symmet- 
ric top about the symmetry axis, 
rad 

solid angle, 4 tt, rad 
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APPENDIX B 


DERIVATIONS OF THE AVERAGE CROSS SECTIONS FOR THE LINEAR MOLECULE 
WITH ZERO AND NONZERO ELECTRONIC POLARIZABILITY 
Zero Electronic Polarizability 


At high J the induced polarizability from equation (18) reduces to 


a iR = 


= — f - — 


(Bl) 


4J^E 


R 


TTi fu 

The 1 partial cross section for a molecule in the Mj state is 


°iR = * P J 


2o iR (Mj)e^ 


1/2 


(B2) 


The average cross section is the properly normalized sum of the over the range of 

Mj values which give positive contributions. This sum can be replaced by an inte- 

gral for high J since the rotational levels are effectively continuous. The normalization 

o 

is 2/(2J + 1) since depends on Mj; this reduces to 1/J for high J. The average 

cross section is then 



Integration with insertion of the limits gives 


_ _ L2 e 2 \ 1/2 jVe - In (Vi + V2) 

* iR 1 7e 


(B3) 


(B4) 
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Nonzero Electronic Polarizability 

At high J the total polarizability is from equation (18b) 


a tR = % + 


4J 2 E 


(m 2 - j2 


R 


(B5) 


The t** 1 partial cross section in the form of equation (19) is 


a tR “ 


it ] 2 e 


J e 


%+ -L( 3 m2-j2) 


4e r j 


a/2 


(B6) 


The average cross section is the integral over contributions corresponding to positive 
that is, 



or 


M q = 0 for a Q > 



36 



The average cross section when the exponent of E-^ is a is 


/ 22 

?tR= *fa) 


Vc + 2 A C - 1 ln ( Vc + 2 +Vi 


2V3 


Vc^~ 


When the exponent is b, the cross section is 


* /3 p 2 e 2 

tR \ 2E R e , 


, 1/2 


Vc + 2 C - 1 ln / Vc + 2 +V3 


2V3 


VT 


(B8a) 


(B8b) 
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APPENDIX C 


COORDINATE SYSTEM FOR MODEL I 

The coordinate system used to describe the interaction between the ion and the linear 
polar molecule is shown in figure 9 (p. 23). The polar angle for translational motion 0 
defines the orientation of the R vector relative to the z-axis. It is measured from +z 
to -z as indicated and goes from 0 to it. The azimuthal angle for translation <£> defines 
the orientation of the projection of R in the x-y plane with respect to the x-axis. This 
angle is measured counterclockwise and goes from 0 to 2tt. The polar angle for rotation 
£ 1 defines the orientation of the negative half of the rod with respect to the z-axis. It is 
measured from +z and ranges from 0 to n . The azimuthal angle for rotation 77 ^ de- 
fines the orientation of the projection of the negative end of the dipole with respect to the 
x-axis. It is measured in the x-y plane in the counterclockwise direction and ranges 
from 0 to 27T. 
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APPENDIX D 


CHOICE OF CONSTANTS AND TRANSLATIONAL ENERGY APPROXIMATION 

The dipole moment \i is presumed to have a constant value along the internuclear 

axis equal to a typical value of 1. 67 Debye units. The representative value chosen for the 

3 

electronic polarizability a Q was 1. 67 angstroms . The moment of inertia Ig employed 
in model I was 3. 85x10"^ gram- centimeter^. This value is 10 times the principal axis 
moment of inertia for the methyl alcohol molecule CHgOH (ref. 10). 

For convenience, all of the relative translational velocity is assigned to the ion. 

This is a good approximation down to where the ion has thermal energy. The ion energy 
e l = ^i v i/2 is simply referred to the center mass by replacing Mj by the effective 
mass of the pair M. The relative energy e is merely Mvj/ 2 . Many IM systems con- 
sist of equal mass molecules and molecular ions where M = Mj/2 = M /2, the "parent" 
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ion case (refs. 6 and 7). The value of M chosen was 5. 12x10 gram for models I 
and n. 



APPENDIX E 


INITIAL CONDITIONS AND EQUATIONS OF MOTION 

The initial conditions for variables and their time derivatives were chosen so as to 
ensure negligible IM interaction at t = 0. This had to be done without entailing prohibi- 
tive running time on the computer. The initial value of R for all cases was 50 ang- 
stroms. Initial values of the four angles 0, and r] were varied as dictated by 

randomizing the rotational plane and rod orientation. The values of the time derivatives 
were judiciously chosen so as to restrict the translational and rotational energies to val- 
ues of interest. 

The Lagrangian for the ion-linear molecule system is 


J 1 = T 1 " V 1 = ~ + r2 ® 2 + r2 sin20 # 


i 2 ) + I (^l +v\ sin 2 ^) + ^° e 


2R 


+ ^2 { C ° S ~ ~ sin ^1 sin 0 ~ cos(?7j - 4>)j| (El) 

R 


The equations of motion follow from 


_d 

dt 




(E2) 


where q = Rj, 0^, and The programed equations in the form to be solved 

for the accelerations are 


" *2 „2 . 2 ^ ’2 
R = R0 + R sin © 



•• • 2 

0 = sin 0 cos 0 $ 



|cos(£j - 0) - jl - cos - <&)J sin sin 0^ 


jcos 0 sin |l - cos (r/j - <£)J - sin^ - 0)| 


MR 


2d! e 2 
o 

MR 5 

(E3) 

(E4) 
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$ = 


(jue/R 2 ) [sin 0 sin £ j sin(rj ^ - $)J - 2MR<f> sin 0(R0 cos 0 + R sin 0) 


MR 2 sin 2 0 


(E5) 


v 1 = . 



sin 0(2<f>R0 cos 0 + R$ 


sin 0 + 2R<& sin 0) + 2Ig? ] ^ sin £ ^ cos £ J 

I B smh j 


(E 6 ) 


^1 = i ? 2 sin cos £ j - ^ e - jsin(£j - 0 ) + [l - cos (77 j - $)] • sin 0 cos 1 -J 


(E7a) 


i b r 


An alternate expression for the 77 ^ equation which was used to check equation (E 6 ) is 


sinfa, - $)sin 0 

77 1 = - afjili ctn . - A — (E7b) 

1 111 ^ sin ?1 

Because of singularities in the coordinate system, an alternate system of equation 
was written in Cartesian coordinates. The problem was transferred to this second sys- 
tem whenever the singularity was approached in the first system and vice-versa. 

The values of constants for which there is meaningful IM interaction can be esti- 
mated by nondimensionalization of equation (El). Let a reference time be 

1 1 2 / juem and a reference distance R R gp be V 21 /M. Then the nondimensional time 
and distance are 


t R~ 


l REF 


(ft = 


R 

R REF 


The Lagrangian becomes 
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L x = + (ft 2 © 2 + <ft 2 sin 2 © 6 2 )+ 1 (i 2 + 77 ! sin 2 ^) 


j (cos(£j - 0) - sin sin 0[l - cos^ - ^ Meo! o 


,4 8/j.I 


The first three parenthetical terms in braces will be of the same order for R equal 
to for times of the order of t = t-^pp. The l/fa^ term will be significant under 

such conditions only if the term Meo; o y4l^8juI is =1. For the chosen values of the mo- 
lecular constants, Rppp is roughly 5 angstroms, that is, several molecular diameters. 
For the same constants, tppp is approximately 2. 5xl0 - * 2 second, the time required 
for an ion traveling 10^ centimeters per second to traverse 2. 5 angstroms. For the hy- 
pothetical molecule, Mea 0 /8jU.I is about 0. 04 so the R"^ term is equal to unity for 
R = 2. 5 angstroms. 
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APPENDIX F 


COORDINATE SYSTEM, DEGREE OF FREEDOM CONSIDERATIONS, AND EQUATIONS 
OF MOTION FOR THE ION-SYMMETRIC TOP INTERACTION 

In the ion- symmetric top coordinate system the radius vector R and angles 0 and 
$ are defined just as in model I. The angle 772 is actually measured counterclockwise 
from the x-axis to the line of nodes (ref. 14), but can be defined as in figure 10. Al- 
though the polar angle ij is identical to it was denoted £ 2 to avoid confusion. The 
additional angle i// 2 defines the orientation of the molecular group CHg with respect to 
the symmetry axis. The angles ?j 2 , l 2 > and ^2 3X6 Ru * er angles, equivalent to the 
angles 6, cp, and t/' of reference 14. 

The ion-dipole potential is the same as in the linear molecule case since the dipole 
is assumed to be "frozen” on the symmetry axis of the top. The angle i// 2 introduces a 
second moment of inertia 1^ about this symmetry axis. The coordinate system is shown 
for a sample prolate top CHgCN. The moment of inertia ratio I^/Ig chosen for the hy- 
pothetical molecule was 0. 2 (prolate top). The angle i// 2 was varied randomly for a 
given orientation of the principal (B) axis as specified by ^ and l 2> 

The value of was varied with 77 2 , ^2’ anc * ^2 so as to insure equipartition 

of energy in the molecule. Initial conditions were assigned so that one-third of the rota- 
tional energy would be present in each degree of freedom. 

The additional degree of freedom introduces the only new energy term in the 
Lagrangian expression. The Lagrangian is 

L 2 = T 2 - V 2 = ~ ( r2 + r2 ® 2 + r2 sin2<E> ^ 2 ) + y (ji\ + ^2 sin2 ^) 

+ — (if; 2 + V 2 cos ? 2 ) 2 
2 

e 2 

+ ii® |cos(| 2 - 0) - sin £ 2 sin 0 [l - sin(? 7 2 - <E>)j| + — - — (FI) 

R 2 L J 2R 4 

The equations of motion for model n are 
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2 a e 

R = R0 2 + R sin 2 0 $ 2 - • |cos(£ 2 - 0) - [l - sin(i 7 2 - <£)] • sin l 2 sin - — — 

MR 3 1 J MR 5 


pe 

2 


0 = sin 0 cos 0$ 2 - - — 


|[l - sin (r] 2 - $)]cos 0 sin l 2 - sin(£ 2 - 0)| 


(F2) 


(F3) 


R 


MR 


i^sin0 sin £ 2 cos (??2 - ®) + 2MR<$ sin 0(R0 cos 0 + R sin 0) 


MR 2 sin 2 0 


(F4) 


^2 


i£e 

R 2 


[sin 0 sin l 2 cosfe, - •)] - 2I B V 2 sin 5 2 cos J 2 * + ^ ^ ^ 


i B sin l 2 


I B sin 


(F5) 


in 


’2 cos h ' I s 


- * (^2 sin h + % sin h cos *2) 

% 


(F6) 


^ 2 - ^2^2 ^2 - ^2 ^2 (F7) 

To avoid singularity difficulties for the symmetric top, an alternate set of equations 
was written in a second Euler system. These equations were utilized in the same man- 
ner as the Cartesian equations written for model I. 
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